Finite-Size Effects in the Field Theory Above the Upper Critical Dimension 
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We demonstrate that the standard 0(n) symmetric </5^ field theory does not correctly describe 
the leading finite-size effects near the critical point of spin systems on a d-dimensional lattice with 
d > 4. We show that these finite-size effects require a description in terms of a lattice Hamiltonian. 
For n — > oo and n — 1 explicit results are given for the susceptibility and for the Binder cumulant. 
They imply that recent analyses of Monte-Carlo results for the five-dimensional Ising model are not 
conclusive. 
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The effect of a finite geometry on systems near phase 
transitions is of basic interest to statistical physics and 
elementary particle physics. In both areas the Hamil- 
tonian 



H = d'^x 
Jv 



(1) 



for an n-component field (p{'x.) in a finite volume V 
plays a fundamental role |]l|. For simplicity we consider 
a d-dimensional cube, V — L"^, with periodic bound- 
ary conditions, (p{x) — L^'^J^i^'Pi^^^^^ ■ The sum- 
mation runs over discrete k vectors with components 
kj — 2'Kmj / L,mj = 0, ±1, ±2, j = \,2,...d, in the 
range —A <kj<K with a finite cutoff A. 

It is generally believed that the leading finite-size ef- 
fects near the critical point of d-dimensional systems can 
be described by H both for d < du and for d > du where 
dti = 4 is the upper critical dimension. Since for d > 4 
the bulk critical behavior is mean-field like, it is plausible 
that the leading finite-size effects for d > 4 appear to be 
describablc in terms of a simplified Hamiltonian 



and imply that recent analyses of the MC data [4-7] in 
terms of the continuum ip"^ theory are not conclusive. 

We shall prove our claims first in the large-n limit 
where a saddle point approach |jl| can be employed. Our 
proof is not based on the renormalization group. We have 
extended the saddle point approach to the finite system 
to derive the order-parameter correlation function 



X 



d'^x < (^(x)(p(0) > 



(3) 



with the statistical weight exp(— TJ). In the limit n — > oo 
at fixed uqu we have found the exact result 



X 



(4) 



We shall denote the bulk critical temperature by Tc. For 
T > Tc, X can be interpreted as the susceptibility (per 
component) of the finite system. In the bulk limit the 
standard equation ||^ for the bulk susceptibility Xb for 
T > Tc is recovered from Eq. (0) as 



rl 



(2) 



involving only the homogeneous fluctuations of the low- 
est (k = 0) mode cpo = L'^^, = L'"^ Jy d'''xip{x). Based 
on the statistical weight exp[—Ho{^)], universal results 
have been predicted B for systems above du ■ For the case 
n = 1, the lowest-mode predictions have been compared 
with Monte-Carlo (MC) data for the five-dimensional 
Ising model [4-7]. Although disagreements were noted 
and doubts were raised in Ref.4, subsequent analyses [5- 
7] based on the Hamiltonian H appeared to reconcile the 
MC data with the lowest-mode predictions. 

In this Letter we shall demonstrate that the lowest- 
mode approach fails for the Hamiltonian H in Eq.(|^) 
for d > 4 and that the leading finite-size effects of spin 
systems on a d-dimensional lattice with d > 4 are not 
correctly described by H. We show that this defect of 
H is due to the (V';')^ term. These unexpected findings 
shed new light on the role of lattice effects for d > d„ 



Xb 



^ = ro + iuon 



Xb 



(5) 



where stands for {2tt)^'^ J d'^k with a finite cutoff 
\kj\ < A. It is convenient to rewrite Eq.(^ in terms of 
^0 ~ foc — Clot where tqc = — 4uon is the bulk criti- 

cal value of To as determined from Eq.(||) ( with Xb^ = 0)i 
and t = {T — Tc) /Tc- Furthermore it is important to scr 
arate the k = part AuquL^'^x from the sum in Eq.(|^ 
After a simple rearrangement we obtain 



_i _ &ro + V(fro)2 + 16uon£^'^(l + S) 
^ ^ 2(1 + 5) 

5ra = aot - A, 

5 = 4^.o«i-'^5][k2(x-i+k2)]-i, 
A = 4Mon [ / k-2 - L-^' k- 



(6) 

(7) 

(8) 

(9) 
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These equations are the starting point of our analysis. 
They are exact in the Umit n — > oo at fixed uqU and are 
vahd, at finite cutoff A, for d > 2, for arbitrary L and for 
arbitrary tq. They are written in a form that separates 
the k = contribution IGuquL^'^ from the effect of the 
k 7^ modes. The latter is contained in S and A. 

In addition to the finite-size effect of the k = 
mode, the k 7^ modes cause two different finite-size 
effects: (i) a finite renormalization of the coupling uon 
due to S which for d > 4 attains the finite bulk value 
Sb = iuon J^ik'^iXb^ + and (ii) a shift of the 

temperature scale due to A which vanishes in the bulk 
limit. These two kinds of finite-size effects were also iden- 
tified by Brezin and Zinn- Justin who argued that for 
d > 4 these effects do not change the leading L depen- 
dence obtained within the lowest-mode approximation. 
These arguments do not depend on n and, if correct, 
should remain valid also in the large-n limit. 

The finite-size effect (ii) comes from A which, for d > 2 
and finite A, has the nontrivial large-i behavior 



0.26 



A AuquA 



d-2 



ai{d){AL)-^ + a2{d){AL) 



2-d 



(10) 



apart from more rapidly vanishing terms. For the coeffi- 
cients ai{d) > we have found 



ai{d) 
02 (d) 



3(27r)'i-2 
-1 



dy[i 
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d/2 



(11) 

(12) 



as confirmed in Fig.l by numerical evaluation of Eq.(|^) 
for d — 3,4,5. Thus, for d > 4, A vanishes as L~^, and 
not as L^-'^ [2,5-7] or as L'"^/^ 0]. This implies that in 
Eq.(|) the zero-mode term proportional to L does no 
longer constitute the dominant finite-size term. 

Our claims are most convincingly examined at bulk Tc. 
Then Eq.(||) is reduced to 



-A + ^/S? -f \&UQnL-'^{l + Sc 
2{l + Sc) 



(13) 



where Sc is given by the r.h.s. of Eq.(||) with re- 
placed by xZ^ ■ We see that the large-L behavior is sig- 
nificantly affected by the A^ term. For large L and d > 4 
we obtain from Eqs.(|l|)and (1^) 



Xc 



L^A 

Auqu 



ai{d)K''-'^L 



(14) 



By contrast, the lowest- mode approximation with A = 
and Sc — Q yields xoc = {Auon)^'^/'^L''-/'^ . This proves 
that the lowest-mode approach fails in the present case. 
We note that the arguments in Ref.2 regarding the finite- 
size effect (ii) are not compelling since they are focused 
on the contributions of individual terms at lowest non- 
zero k rather than on an analysis of the summed effect of 
these contributions. 
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FIG. 1. AL-dependence of Ao = A/(4MonA''"^) with A 
from Eq.(^) for d = 3, 4, 5 (solid curves). The dashed lines 
represent Eq.(0 with ai(3) = 0.27706, ai(4) = 0.08333, 
ai(5) = 0.02443, and 02 (3) = 0.22578, 02(4) = 0.14046, 
02(5) = 0.10712. The arrows indicate the large-AL limits. 

Furthermore we see that the L'^~'^ power law in Eq.(p^) 
differs from the L'^^'^ power law obtained from the exact 
solution of the n- vector model on a lattice for n — > cxd |^ 
and of the mean spherical model on a lattice [|o| . This 
proves that the field-theoretic Hamiltonian H in Eq.(j^) 
does not correctly describe the leading finite-size effects 
of spin models on a d-dimensional lattice with d > 4, at 
least in the large-n limit. 

In the following we show that this defect is due to the 
{S/'^Y or k^vsktjf-k term of 



+UqL ^'^ ^ ('^k'^k')('<5k"'^-k-k'-k") (15) 



kk'k" 



with iy9k = jyd!^^^ *'''^iy9(x). Instead we consider a lat- 
tice Hamiltonian H{ipi) for n-component vectors (pi with 
components ifia, — 00 < ipia < 00, a = 1, n, on lattice 
points Xi of a simple cubic lattice with volume V = L'^ 
and with periodic boundary conditions. We assume 




(16) 



where is a pair interaction and d is the lattice spacing. 
In terms of (p\^ — af^ e~^^'^^ipj the Hamiltonian II has 
the same form as Eq.(n5[) but with ro -I- k^ replaced by 
fo + 2(5J(k) where (5J(k) = J(0) - J(k) and 
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J(k) = (a/L)' 



-'ik-(xi— Xj) 



(17) 



The k values are restricted by —it /a < kj < tt/S,. In 
the large-n hmit at fixed iion the susceptibihty x — 
^{a/LYY.^,J < V^V3 > is determined by Eqs.(|)-(|) 
with replaced by 2(5 J(k). The large-L behavior of 
the crucial quantity A is for d > 2 

A = 2uon(^ J [5 J(k)]"i - L-'' ^[W(k)]-i) (18) 



k#0 



-1„ ^J\r2-d 



47/071 Jq a2{d)L 



(19) 



which for d > 4 differs from that of the continuum ver- 
sion A, Eq.(p^, where 2(5 J(k) was approximated by k^. 
This approximation turns out to be the unjustified for 
(i > 4. Eq.(p^) is valid for short-range interactions where 



Jo = ^(a/L)'^^(j.,/a')(x. 



(20) 



is finite. As a consequence of Eq.(|19|), the leading L de- 
pendence of X Tc is for d > A 



Xc-i(uon)-V2(i + 5,^)V2i<i/2. 



(21) 



This agrees with the i'*/^ power law of the exact solu- 
tion of the lattice models of Refs. 9 and 10 for d > A and 
with the lowest-mode result xoc = ^{uon)~^^^L'^^^. We 

see that at Tc for d > A the k 7^ modes of H do not 
change the leading exponent d/2 of the lowest-mode ap- 
proximation. Nevertheless they produce a finite change 
of the amplitude of Xc through = uqu /^[(^Jfk)]^^. 
Furthermore we conclude from Eqs.(|lJ) and (|2^) that 
for d > A the lattice Hamiltonian H yields significantly 
different finite-size effects compared to those of H . 

An analysis of the temperature dependence of x(i, L), 
Eq.(^, and of x{t,L) shows that for (i > 4 finite-size 
scaling in its usual form is not valid, as expected j9j, but 
we find that it remains valid in a generalized form with 
two reference lengths. The asymptotic scaling structure 
of X is 



x(^,L) = i^/'^4(^(L/eo)'/^(i/^o) 



(22) 



where ^0 is the bulk correlation-length amplitude and 
fo = [4uonJo'^(l + S'^)~^]^/('^~^) is a second reference 
length. The (i-dependent scaling function reads 

P^{x,y) = 2J^^[5{x,y) + ^Mx^WT^Y' (23) 

where 6{x,y) — x — a2{d)y. In the lowest-mode approxi- 
mation the term —a2{d)y is dropped which implies that 
the leading finite-size term for t > becomes incorrect. 
Thus, for t > 0, the lowest-mode approach fails for the 



lattice model (and also for the continuum model whose 
scaling function turns out to be non-universal). 

In the following we extend our analysis to the case 
n — 1 which is of relevance to the interpretation of MC 
data of the five-dimensional Ising model [4-7] . We shall 
examine the susceptibility x ^nd the Binder cumulant U, 



X 



d'^x < V3(x)^(0) >= L'^ < $2 >, 



1 



C/ = 1 - ^ < > / < >^ 



(24) 
(25) 



within the ip^ model, Eq.(|I]), including the effect of the 
k 7^ modes in one-loop order. For (i > 4 at finite cut-off, 
the perturbative finite-size field theory 0,|2| is applica- 
ble near Tc without a renormalization-group treatment. 
The averages are defined as < $™ >= (i<i>$'"P(<i>) 
where P{^) = J T>ae~^ is the order-parameter dis- 
tribution function with (j(x) — <jo(x) — $ representing the 
inhomogeneous fluctuations |1^. From Refs. 11 and 12 
we derive the L dependence at Tc in one-loop order 



Uc 



where 
Vo 



r^ffTd/2 eff-1/2 

'0 - roc + 12uoSi{roL) 



lAAulMlS2{r^L), 



tql = roc + 12uoMo , 

Jg°° (iss™ exp(-iys^ - s^) 



Jp°° ds exp(- 



1 ? 



(26) 
(27) 

(28) 
(29) 
(30) 
(31) 
(32) 

(33) 



In this order the critical value rpc < is determined im- 
plicitly by the bulk limit {rl'^^ = 0) of Eq.(|29l), 



12wo/i(-2roc) + 36Moroc/2(-2roc), (34) 
k^)^™. The finite-size effect of the 



where /^(r) - 
k 7^ modes enters through 



(r + k^ 



(35) 



For large L we have tol = — 2roc + 0{L and 

r'o" = - 12uo[hi-2roc) - 5i(-2roe)] 

+ 36uorocM~2roc) - ^2(-2roc)] + 0(1-''). (36) 

Similar to A in Eqs.(^ and (|lO|), the parameter rg^-^ < 
vanishes as (rather than as L'^^'^) for d > 4, thus 
y^^^ < diverges as L^'^^'^^/'^ (rather than vanishes as 

^(4-rf)/2^ for d > 4. Since d2{y) 2//4 and ^^{y) - 

2/^/16 for large negative y this implies that Xc di- 
verges as U^^^ and Uc attains the large-L limit 2/3. We 



3 



conclude that the L''/^ power law for xoc and the value 
Uoc = l- ^^i{0)/MO? = 0-2705 predicted |] for n = 1 
within the lowest-mode approach are incorrect. From 
Refs. 2 and 12 we infer that analogous conclusions hold 
for general n > 1. 

These unexpected results show that the widely ac- 
cepted arguments in support of the asymptotic correct- 
ness of the lowest-mode approximation above the upper 
critical dimension in statics [1,2,5-7,13-19] and dynamics 
[1,20-23] are not valid and that recent interpretations [4- 
7] of the Monte-Carlo data of the five-dimensional Ising 
model in terms of predictions based on the Hamiltonian 
H, Eq.(^), are not conclusive, in spite of the apparent 
agreement found in Refs. 5-7. 

Guided by our exact results in the large- rt limit, we 
propose a solution to this puzzle by replacing the field- 
theoretic (^^ Hamiltonian H, Eq.([^), by the lattice (^^ 
Hamiltonian H, Eg. (p^) , with n = 1 for the compari- 
son with the five-dimensional Ising model. This involves 
a reexamination of r^^^ , Eq.(|3^), with replaced by 
2[J(0) — J(k)]. We anticipate that the resulting value for 
the Binder cumulant Uc of the lattice model will be close 
to (or possibly identical with) that of the lowest-mode 
approach. This expectation is based on our result (for 
n — > oo) that at Tc the lowest-mode approach yields the 
correct leading finite-size exponent of Xc- In addition, 
however, a detailed analysis of non-asymptotic (finite-L) 
correction terms is required which, for the Hamiltonian 
i/, are expected to be different from those employed in a 
recent analysis based on H [^. 

We summarize our findings for the continuum and lat- 
tice versions of the (p'^ model for d > 4 as follows: Lat- 
tice effects manifest themselves not only in changes of 
nonuniversal amplitudes but also in changes of the expo- 
nents of the leading finite-size terms as compared to the 
exponents of the continuum (p'^ model. The lowest-mode 
approach fails for the continuum ip'^ model, and also for 
the lattice model for t > 0. Therefore the values for the 
amplitude ratios derived previously ^ cannot be justi- 
fied on the basis of the (p'^ continuum theory. For the 
lattice Hamiltonian, however, the lowest-mode approach 
is qualitatively justified at Tc for n — > oo, at least for Xc- 
We conjecture that this is the reason for a fortuitous (ap- 
proximate) agreement found between MC (lattice) data 
[4-7] and the lowest-mode predictions 0]. Further work 



is necessary in terms of the ip'^ lattice model, Eq.(16), to 
fully establish our conjecture. 

We also anticipate lattice and cutoff effects on lead- 
ing finite-size terms at d = du- This is relevant to future 
studies of tricritical phenomena at d = 3, e.g., in '^He-^He 
mixtures [Q and to MC simulations for lattice models 
of elementary particle physics at d = 4. 
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